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Abstract 

We discuss Hall effect and power dissipation in chiral p-wave superconductors near Kosterlitz- 
Thouless transition in the absence of applied magnetic field. In bound pair dynamics picture, 
nonzero Hall conductivity emerges when vortex-antivortex bound pair polarization has a compo- 
nent transverse to the direction of external perturbation. Such effect arises from the broken time 
reversal symmetry nature of a chiral p-wave superconducting state and does not require an ap- 
plied magnetic field. A frequency-dependent matrix dielectric function e(u) is derived to describe 
the screening effect due to the pair polarization. Quantities related to the Hall conductivity and 
power dissipation, denoted as ej 1 and 9(— e^ 1 ), are investigated in frequency and temperature 
domain. The imaginary part of the former can show peak structure and sign reversal as a func- 
tion of frequency close to transition temperature, as well as in the temperature domain at various 
fixed frequencies. The latter shows peak structure near transition temperature. These features 
are attributed to pair-size-dependent longitudinal and transverse response function of bound pairs. 
Consequences due to free vortex dynamics and the resulting total conductivity tensor a are also 
discussed. 

PACS numbers: 73.50.Jt,47.32.C-,74.20.Rp 
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I. INTRODUCTION 



Zero-field Hall effect in chiral p-wave superconductors (SCs) has drawn much attention in 
literature recently.- - - Because of the nature of broken time reversal (T) symmetry, a nonzero 
Hall conductivity can be possible in a chiral p-wave SC. Indeed, it has already been shown 
that spontaneous Hall effect could arise from the intrinsic angular momentum of Cooper 
pairs^ as well as from the spontaneous surface current.- More recently, Hall conductivity 
due to impurity effect^ or to multiband SC structure^ was also studied, which could 
give possible explanation to the observed polar Kerr effect in the superconducting state of 
Sr 2 Ru0 4 .^ 

In this work, we address the zero-field Hall effect in a chiral p-wave SC originating from 
another mechanism, namely the vortex dynamics near Kosterlitz-Thouless (KT) transition. 
In two-dimensional (2D) superfluid (SF) or SC films, quantized vortices are realized as 
topological defects in the condensates, whose dynamics has been one of the key ingredients in 
understanding 2D phase transition phenomena.- A few decades ago, Kosterlitz and Thouless^ 
suggested a static theory to relate a phase transition observed in superfluid 4 He filrn^ 
to vortex- antivortex pair unbinding process across a transition temperature Tkt- In this 
picture, the logarithmic vortex-antivortex interaction is screened by smaller pairs and is 
renormalized to 1/eof its bare value Kq for temperature T < Tkt- The length-dependent 
dielectric constant <T is used to describe the static screening of pair interaction. When T > 
Tkt, there exists a finite pair size £ + such that the interaction becomes vanishingly small. 
Consequently, the pair unbinds and free vortices emerge; superfluidity is then destroyed. 
Soon after that, Ambegaokar, Halperin, Nelson, and Siggia (AHNS) 1 ^— combined this static 
theory with Hall and Vinen's dynamical description of vortex motion^ 4 - to give an analysis 
of the dynamical effect on the phase transition. Concisely speaking, the renormalization 
process in the static theory^ is truncated by vortex dynamics with a characteristic length 



7\d(o;) = a/14-D/oj instead of going to its completion.— Here D is the diffusivity of vortex 
movement and u is the driving frequency. This results in broadening transition observed in 
4 He SF films 1 ^ as well as in charged Fermi systems such as high-temperature SCs.— 

It is instructive to explore any physical consequence stemming from this KT transition in 
a broken T symmetry state. Possible experimental candidates of broken T symmetry state 



could be superconducting Sr 2 Ru04 (Ref. Il8l ) or He- A phase thin film,— in which pairing of 



chiral p-wave type is expected. Indeed, in literature some theoretical works have been done 
to investigate new features specific to SCs with pairing of this type near KT transition.— 

In this work, we consider a 2D p-wave pairing state with d- vector d = z(p x ± ip y )/pF 
where z is the unit vector normal to film surface, p x and p y denote the x and y component 
of the relative momentum p of a Cooper pair, and pp is the Fermi momentum. Assuming 
isotropic Fermi surface, two kinds of pairing fields can be obtained: A( + ~)(p, p)/A QO = 
rt + -\p)e- i «( Px + ip y )/ PF + ft~\p)e i ' t> ( Px -ip y )/p F and A^\r,p)/A OQ = f[ ++ \ p ) e ^(p x + 
ip y )/pF + f2 ++ \p) e3l( ^(Px — ip y )/ P F with asymptotic behavior at large p being f[ + — > 1, 
f2 + ~' ) — > 0, fx' — > 1, and / 2 (++ ^ — 0.— ~— Here p = p(cos0, sin0) is the spatial coordinate 
and Aqo is the modulus of the energy gap in the bulk. From these pairing fields, we can 
identify two types of integer vortices called (H — )-vortex and (++)-vortex respectively. 

Because of spontaneously broken T symmetry, these two types of vortices are not 
equivalent.— ~— In particular, their Hall and Vinen coefficients^ do not share the same 
value, i.e., B^ ++ ^ ^ B^ ) and B'^^ ^ B'^ + ^ (see section II A). This results in a nonzero 
"convective" term in a vortex pair polarization Fokker-Planck equation in addition to the 
conventional diffusive terms, while in its s-wave counterpart such convective motion does 
not enter the dynamics.— The relative strength of convection is quantified by a convective 
ratio xq in this paper. It is due to such distinct feature that pair polarization transverse 
to the driving force field becomes possible even without applied magnetic field. A nonzero 
vortex-dynamics-induced Hall conductivity cr± then follows naturally. 

The main result of this work is that in the bound pair dynamics description, we obtain 
nonvanishing Hall conductivity a± and ac conductivity an near the KT transition. One of 
the interesting features in the Hall conductivity is that strong positive peak and sign changes 
in o;9ft(<7_i_) are observed at suitable frequency region above the transition temperature, as 
well as above Tkt in temperature domain at fixed frequencies. On the other hand, oj§t.(a\\) 
is shown to have similar features as in AHNS's results. We note that the shapes of two 
length-dependent response functions Q\\(r,u) and Q±(r,u), which corresponds respectively 
to the longitudinal and transverse response of bound pairs with separation r to external 
perturbation with frequency u, play a determining role on the behavior of a± and cry. We 
also discuss the contribution of free vortex motion and the resulting total conductivity tensor. 

The paper is organized as follows: In section II, we generalize AHNS's vortex dynam- 
ics in the chiral p-wave context. To describe the vortex-antivortex bound pair dynamics, 

3 



the above-mentioned response functions G\\(r,w) and Q±(r,u) are derived from the Fokker- 
Planck equation governing the pair motion. Together with the free vortex contribution, we 
arrive at a matrix dielectric function e(u) which describes the total screening effect under 
time- dependent perturbation. In section III, we investigate the frequency and temperature 
dependence of the conductivity tensor a(u>) constructed from e(u), treating the bound pair 
and the free vortex contribution separately. The behavior of total conductivity a is also 
discussed. A summary and remark are given in section IV. Finally, we discuss analytic ex- 
pression of e^ 1 ), ^(ej 1 ), and ej 1 ) in opposite limit of the convective ratio xq <C 1 
and x ^> 1 in Appendix A. 

II. MATRIX DIELECTRIC FUNCTION 
A. Equations of motion for single vortex 

To construct a matrix dielectric function, we consider a neutral SF film system resembling 
that employed in AHNS's dynamical theory, with film thickness of order the superconducting 
coherence length and linear dimension L (W) along x (y) direction. L is very large and W 
is large but finite. Vortex core motion relative to the local superfluid velocity leads to a 
Magnus force -Fjyn^ 

F»=^k(„? -„«). (1) 

In the above equation, and vj? are the velocity of the i-th vortex core and the local 
superfluid flow excluding the diverging self-field of the i-th vortex respectively. rij(= ±1) is 
the vorticity of the i-th vortex. p° is the bare areal superfluid mass density, which is defined 
as the three-dimensional superfluid density integrated across the film thickness, m* is the 
mass of the constituting particle, which is equal to the mass of a Cooper pair. 

The SF film is driven by a vibrating substrate. A vortex core moving relatively to the 
substrate experiences a vorticity-dependent drag force F^— 

F« = BW (v n - vg>) + B'^z x (v n - v^) , (2) 

where v n is the moving substrate velocity. and B'^ denote the vorticity-dependent 

drag coefficients originating from interactions with the substrate and with thermally excited 
quasiparticles and collective modes. Here B^ (B'^) = Z?( ++ ' (B'^ ++ ^ ) for rij = 1 or B^ + ^ 
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for Hi = _i. These quantities have been obtained for a three-dimensional clean 
s-wave SF/SC with isotropic Fermi surface,— and are related to the relaxation time r for 
the Caroli-deGennes-Matricon mode^ in the vortex core. Due to broken T symmetry, the 
relaxation time for the mode in the (++)-vortex and (H — )-vortex core, and thus the 
values of their drag coefficients, are different, resulting in the vorticity-dependent drag force 
on the vortex core. These drag coefficients for 2D clean SF/SCs with cylindrical Fermi 
surface can be inferred from the three-dimensional result easily. 

When two forces balance F^ + F$ = 0, the z-th vortex velocity is expressed by 



vf = m2TrKoDW-jrz x (v n - w») + C®(v n - v^p) + v® + r/«(t), (3) 



where 



(2hirpV/m* - B'®) 2 + 



= ksT — - ; ^ ; — ~ 2 , (4) 



C (i) = j _ 2h *Ps 2h7rp° s /m* - B'® 

m* {2tmp°Jm* - B' {i) ) 2 + B^ 2 ' 

Here K = h 2 p°J (m* 2 k B T) with k B being the Boltzmann constant. t]W(i) are fluctuating 
Gaussian noise sources incorporated to bring the vortices to equilibrium. Their components 

13, the local superfluid velocity v!p is 



satisfy {r}%> (t)r)X> (V )) = 2D®8 ij 8 llv 8{t-t'). From Ref . 
related to the spatial average superfluid velocity u s and the positions of individual vortices 
r,j by Vs = u s + (h/m*)zx^2j J L i 7ijV(3(7*i, rj). Here C5(r, rj) is the Green function satisfying 
V 2 C5(r, rj) = 2irA(r — rj) and boundary condition 0(r, rj) — on the edges. The function 
A(r — rj) is localized in a region around rj with radius of the order of coherence length (we 
could say that the function A(r — rj) is a delta function in the coarse-grained scale). Far 
away from the edges, (3(r,rj) ~ ln(|r — rj\/a,o). The spatial average of V(3 turns out to 
be zero and thus u s represents the spatial average of v$ . Finally, the time evolution of u s 
obeys dtu s = — £\ rii(27rh)/(LWm*)z x w£ , reflecting the fact that the average superfluid 
velocity in the x direction changes by 2nh/ (Lm*) when a vortex with = 1 moves across 
a strip with width Sy. 

In a charged system, if we follow Kopnin's description,— the driving force on the vortex 
due to a transport current j tI is a Lorentz force F-jp = nj$o[jtr x z sgn(e)]/c where $o — 
hmcj\e\ and e is the electron charge. This force is balanced by the force from environment 
■Fenl = —rj^v^ + riirf^'z x v{/ . If we set v n = in Eqs. (J2J) and (j3J), the results derived in 



a neutral system can be carried over to a charged system by the translation 77 M and 



B. Fokker-Planck equation for bound pair contribution 

We consider the polarization of a test vortex-antivortex pair whose constituting vortices 
interact via a screened logarithmic interaction. The pair is under the influence of an in- 
finitesimal oscillating external field 5E. The Langevin equation for their relative coordinate 
r can be obtained by subtracting Eq. ([3]) from each other for opposite vorticity 



Here D = (£>(++) + £>(+-))/2 and C = (C(++) - C^ + ' ) )h/(2m*). U{r) is one-half the 
dimensionless potential energy of the pair and the Gaussian noise now satisfies (77^77^) = 
AD8n V 8{t — £'). The potential energy is given by 



In the above equation, the first term on the right hand side describes the logarithmic interac- 
tion screened by the Kosterlitz dielectric constant e(r). 2/jo i n the second term is related to 
the energy required to create a pair with separation ao- In the last term, SE has dimension 
of velocity and acts as the perturbation. In the integration limit ao is a length scale related 
to the size of a vortex core, and r is the pair separation. 

We can see in Eq. fl6]) that in addition to the conventional diffusive terms depending 
onD, a convective term proportional to C also enters the dynamical equation. While the 
strength of the former is proportional to the average of that of the latter is related to 
the difference of between opposite vorticity. We emphasize here that such a convective 
pair motion is one of the special features for a system with unequal opposite vortices, and is 
thus absent in an s-wave SF/SC since = in that case. Given this nonzero C, the 
pair polarization is tilted away from the direction of the force field — VL r (r), and has both 
longitudinal and transverse components even without applied magnetic field. 

The Fokker-Planck equation corresponding to Eq. ((6]) is given by 



dr 
~dt 



4nDK VU(r) + 2Cz x W(r) + 77. 



(6) 





dT(r,t) 
dt 



AnDKoV ■ (rW) - 2CV • (Tz x W) + 2£>V 2 r, 



(8) 
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where T(r,t) is the density of pairs per unit area of separation. We take the time- 
independent state T to be exp(— 2nK Uo)/aQ where Uq = U\se=o- Now we follow the 
standard procedure,— letting r = T + ST and keeping terms to first order in SE. In 
frequency space, we obtain 



iuST(r,u) = -AirDK V ■ 
+2CV 



777. 

~h SET ° 
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—z x 5ET } 



2DV 2 ST. (9) 



We employ the expansion ST(r, u) = SF m (r, u)e me , where 9 is the angle measured from 
SE/\5E\ to f in anti-clockwise sense. Only ST m with m = 1 and m — — 1 are coupled to the 
external field. 

We define a convective ratio Xq = C/(2irK D), which measures the relative strength 
between convection and diffusion, and introduce an ansatz 5T±i(r, u) = 2ttKq (1 ± ixq) 5ETq 
g±(r,u)rm* /(2h). Together with an approximation^ 2'KK^j'e = 4 (which is valid near the 
transition) and a change of variable z 2 = —iur 2 / {2D), the equations of g±(z) in the ansatz 
corresponding to angular momentum m = ±1 are given by 

= - (z 2 + 4 ± 4ix ) g±(z) - zg' ± (z) + z 2 g" ± {z) + 4. (10) 

We are then able to write down the change of distribution function 5T = 5Tie ld + 5Y^ie~ lB 
in terms of 5E explicitly 

777 

ST(r, to) = 2tiK T — [SE ■ rQ\\ (r, w) - SE ■ (r X z)G ± (r, to)] , (11) 

where 

C?ll(r,w) = ~ [G+(r,u) + G-(r,u)] , G ± (r) = -- [G + (r,u) - G-(r,u)] , (12) 

and G±(r,u) = (1 ± ix )g±(r, uj). In Eq. (ITTj) . Q\\(r,ui) takes the role of pair-size-dependent 
response longitudinal to the driving field SE and Qj_(r, u) is the transverse response func- 
tion. In the limit of vanishing xq, G\\{t,oj) reduces to the s-wave SF/SC result g s (r,u) w 
IAD /{IAD— iur 2 ) while Gj_(r, u) becomes identically zero. In Eq. (1121) . we find that G\\.j_(r, u) 
depend on two quantities G±(r,oj) which describe pair motion with angular momentum 
m — ±1 respectively. 

We describe G\\ t x(r,cu) first. In Fig. [T] we plot G\\(r,u) and G±(r, u) as a function of 
/ = ln(r/a ) for fixed frequency ual/(2D) = 10~ 2 for (a) x = 0.1, (b) x = 0.5, (c) x = 2. 
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FIG. 1. (Color online) The longitudinal and transverse response function G\\(r, uj) (left column) 
and G±(r,oj) (right column) are plotted as a function of I = ln(r/ao) for xq = (a) 0.1, (b) 0.5, and 
(c) 2.0 respectively. The real and imaginary part are indicated by blue dashed lines and red solid 
lines respectively. G\\(r, uj) is qualitatively similar to the s-wave SF/SC result. G±(r, uj) is a new 
feature with a peak in the real part and a dip-recouping shape in the imaginary part around 



Blue dashed lines and red solid lines represent their real and imaginary parts respectively. In 
the left column of Fig. [TJ the longitudinal response function Qn (r, uj) is qualitatively similar 
to the s-wave SF/SC result g s (r,uj) even for finite Xq. It has a step-function-like real part 
and a delta-function-like imaginary part concentrating near (/* marks the position where 
^(^-iJ changes sign; see Appendix A). The response function G±(r, uj) is a new feature in 
this model. In the right column of Fig. [H §t[G±(r, uj)] has a peak structure. This means 
that neither smaller nor larger pair gives response in transverse direction. Only pairs with 
characteristic pair size /* can give rise to the Hall effect. Besides, Q[G±(r, uj)] shows negative 
dip shape for / < /*, and then recoups when I > /*. It turns out in later section that such 
a dip- recouping antisymmetric shape in 3((?_i_) about /* plays a determining role in many 
features of the Hall conductivity. These features become more significant when xq increases. 
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FIG. 2. (Color online) The response function G+(r, w) (left column) and G-(r, uj) (right column) 
corresponding to angular momentum m = ±1 motion are plotted as a function of I = ln(r/ao) for 
xq = (a) 0.1, (b) 0.5, and (c) 2.0 respectively. The real and imaginary part are indicated by blue 
dashed lines and red solid lines respectively. Asymmetry between the two columns becomes more 
prominent when xq increases. 



In Eq. ( fl~2l) . we can see that Q\\(r,iu) depends on the average of G±(r,u), i.e., the average 
response with +1 and —1 angular momentum, and Q±(r, u) is related to the difference 
between them. This means that any asymmetry between m = ±1 motion gives rise to 
nonzero transverse response. In Fig. [2] we plot G±(r,oj) as a function of I using the same 
set of u and Xq as in Fig. [lj Again, blue dashed lines and red solid lines represent their 
real and imaginary parts respectively. We observe that asymmetry between m — ±1 motion 
grows with increasing xq. For small xq in Fig. EJ^a), the motion described by G+(r,oj) and 
G-(r, u) is only slightly asymmetric, and they look like the familiar curve g s (r,u). When 
xq increases in Fig. [2](b) and (c), the asymmetry between the left and right column becomes 
more and more significant. It is worth noting that, when we subtract from each other the 
real part of the functions for opposite m, we can obtain the dip-recouping shape of Q(Q±). 
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Equation (jTOj) can be solved exactly to give the length-dependent response functions, 



12 by 



but here we can employ approximate solution to g±(r,uj) in a manner similar to Ref. 
neglecting g± and g'± in Eq. ( 1TU1) . Then Q^± become 

r ( \ ~ 1 ( l + ^o \-ixq \ 

^ll^'^-oh : — TRK~- V 1 ■ 9 A n — _ ' (13) 

z \ 1 — loor^/ XD + zrco 1 — ioor z / AD — ixq J 

Q±(T,u) ~ 4 ( ^ ^ 1 , (14) 

2 \ 1 — iujr 2 1 XD + ixo 1 — iujr 2 j XD — ixo J 

where A is a factor to fit the exact curve. It is selected to be 14 when xq = and is of order 10 

for a wide range of xq. From Eqs. ( fTBl and (114j) we can see that for a small convective ratio 

1 > xq > 0, there are poles given by 1 — icor 2 / (XD) ± ix = 0. This means that the pair size 

with which a pair gives a strong response differs from the standard result ro = a/14D/o; by a 



length of order ±£oy XD/uj for motion with m = ±1, creating the asymmetry demonstrated 
in Fig. [2l For a large convective ratio Xq > 1, such pole is removed for m — —1. Therefore, 
a broad and flat curve is expected [Fig. [2](c) right column]. 



C. Matrix susceptibility Xb f° r bound pair contribution 

In this subsection, we introduce a susceptibility matrix Xb for the bound pair polarization 

(15) 



Xb 



( xl xi 



\~Xb x b/ 



whose components are defined as 



2ixh f 2 ( r 5V \ ± 27rh f l2 „ ( r 



Xl = — / d 2 r - ~ , Xb=— d 2 rz ■ - X ^- . (16) 
Ab m* 7 V 2 ^/ Ah m* J \2 5E J K ' 

They are so defined that the real parts of xt an< ^ Xb are positive when xq > 0. Using 
Eqs. ( 1TT]) and f fl6|) . we arrive at the expression 

xl = dr-g^r), Xb L = J ^ dr-G^r). (17) 

The infinitesimal external field redistributes the pair polarization by an amount 5P^ accord- 
ing to 5Pt, = Xb$E. The definition in Eq. ( f!6|) means the external field tilts the polarization 
in clockwise direction when xq > 0. It should also be noted that the integration are per- 
formed from a to £ + where £ + is a coherence length with behavior £ + — > oo for T < Tkt 



and £ + ~ exp[2n /(by^T/T KT — 1)] for T > T K t- & is a non-universal positive constant of 
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order unity. From Eq. (jTTj) it is now clear that while the longitudinal polarization Q\\ takes 
the role of the function g s (r,u) in AHNS's theory, the transverse polarization Q± has no 
simple analogy in s-wave SF/SCs and is specific to systems with finite convective ratios. 
Integrating the two pair-size-dependent response functions with the weight function de/dr 
gives the susceptibility matrix elements from the bound pair contribution. 



D. Free vortex contribution and total dielectric function 



As for the free vortex contribution, from Eq. ([3]) we can obtain the equation for polar- 
ization Pf for plasma of free vortices under the spatial averaged driving field (E). The total 
free vortex density is given by rif = N/ (LW) with N being the total number of free vortices 
and LW being the area of the film. In frequency space, we have 



2nh 

— iuPf = rif 

m* 



2nDh p° s m* — „ 

\E) — Cz x (E) 



m* ksT h 
If we define the susceptibility for free vortices as Pf = Xi(E), we have 



(18) 




Xf=^( 
— iuj 1 

where 70 = A7r 2 n{DK . On the other hand, the free vortex density is related to the coherence 
length by rif = 2nC\£ > + 2 where C\ is a positive constant of order unity- 
Together with the bound pair contribution, the total dielectric function can be obtained 
as e(cu) = 1 + Xb + Xf an d its inverse reads 

where e^ 1 = -{ X £ + Xt)/[0- + xl + x\f + {xt + xi?} and e" 1 = (1 + X \ + Xf)/[(1 + xl + 
Xf ) 2 + (Xb + X~t) 2 ]- We are now in position to present our results on the Hall conductivity 
and power dissipation due to this dielectric function. 



III. RESULT AND DISCUSSION IN CHARGED SYSTEM 

It was discussed that the conductivity tensor a in a charged system was related to the in- 
verse dielectric function e -1 discussed above.— We can understand the relation by considering 
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the total current under the influence of a driver coil electric field £.— The vortex-modified 
total current jf tot is related to the field by £ = —iuLkej to t where Lk is the sheet kinetic 
inductance. This relation is generalized to our model with a matrix dielectric function e. If 
we use the sign convention that the normal state electron Hall conductivity with magnetic 
field pointing in z direction is positive, cr{w) and e~ 1 (w) is related by 



In particular, we may write down ojLfc5R(o"j_) = ^(e^ 1 ), wLfcQ^ajJ = 9R(— ej 1 ), and 
u>Lk$l((T\\) = < ^s(—e^ 1 ). Among these expressions, the first two are associated with the 
real and imaginary part of the Hall conductivity, while the last one is related to the power 
dissipation P because P oc 9ft(o"||) oc S(— e^ l )/uj. 

In the following subsections, we investigate the frequency and temperature dependence 
of ^(e^ 1 ), 3?(— e]^ 1 ), and 3(— e^ 1 ) respectively. 

A. Hall conductivity from bound pair contribution 

We consider bound pair contribution to the Hall conductivity first by ignoring the terms 
x|' and xt in e ~ l ■ Plots of the negative real part (blue circles) and imaginary part (red 
squares) of e^ 1 versus (a) logcU and (b) t are presented in Fig. |3j Here, ZJ = coal/(2D) 
is the scaled frequency, and t = T/Tkt — 1 is the reduced temperature. In this and the 
following figures, we use the renormalization group flow equations <9;(2 — X{)~ 1 = 47r 2 y z 2 and 
d x Yi = XiY h H with initial conditions X = 2 - [2 - X^ KT) }/(1 + t) and 2itY = 0.1 related 
to the bare interaction strength K and chemical potential /xq respectively. At transition 
temperature, Aq KT ^ = Xq(T = Tkt) ~ —0.215 which is obtained numerically. Here, it 
suffices to notice that Xi is given by Xi = 2 — irKi where Ki = 'e~ 1 K is the renormalized 
interaction strength, and Yi is related to the the number of bound pairs found with pair size 
r = a^e 1 . Also, we use the convective ratio Xq — 0.1 in the plots for illustration purpose. 

We first focus on the discussion of the frequency dependence of ej 1 at some fixed tem- 
peratures in Fig. Ela). In the low temperature phase t < 0, both 9R(— ej 1 ) and ^(ej 1 ) are 
positive and increase steadily with frequency. When temperature increases to the high tem- 
perature phase t = 0.045, ^(ej 1 ) shows a positive-valued peak followed by sign reversal at 
higher frequency indicated by the red horizontal double-headed bar. Meanwhile, 5i(— ej 1 ) 




(21) 
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FIG. 3. (Color online) Bound pair contribution. The negative real part (blue circles) and imaginary 
part (red squares) of ej 1 as a function of (a) logcJ and (b) t = T/Tkt — 1- The other parameters 
used are xq = 0.1, 2ttYq = 0.1 and Xq ~ —0.215. In (a), above the transition temperature, 
^(e]! 1 ) has peak structures and sign reversals. ej 1 ) has a sharp increase when frequency 
increases. In (b), ^(ej 1 ) shows peak structure and sign reversal above t = and K(— ej 1 ) has a 
peak. The features broaden and move to higher temperature when frequency increases. 

increases sharply around the peak of ^(e^ 1 ). These features move to higher frequency side 
when temperature increases further (t = 0.090). For temperature dependence at fixed fre- 
quencies in Fig. [3](b), the sign reversal and peak structure in ^(ej 1 ) are also observed when 
temperature varies above Tkt- Meanwhile, 5R(— e^ 1 ) shows a simple peak structure across 
Tkt- When frequency increases, the structure broadens and moves to higher temperature 
region. 

Sign anomaly in Hall conductivity has been observed in various superconducting systems 
such as high-temperature SCs,— r — and it is known that single vortex "upstream" motion or 
change of sign of charge carriers could give rise to such phenomenon.— Here in our model 
the sign change in ^(ej 1 ) stems from the vortex-antivortex pair unbinding process. For 
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a vortex-antivortex pair moving "downstream" with same speed, total transverse electric 
field is canceled and no Hall effect can be observed. Now that the constituting vortices are 
not equivalent to each other, they can respond differently under the influence of transport 
current, and a net transverse electric field follows. 



Pair polarization Pair polarization 




Positive Hall signal Nagative Hall signal 



FIG. 4. (Color online) Schematic diagram of vortex-induced electric field E^ under transport 
current j tr with convective ratio xo > 0. The black solid arrow is the transport current j tr . 
The blue and red solid arrow represent the pair polarization of pair with size I < I* and I > Z* 
respectively. The dashed arrows show the electric field generated by pair motion in the respective 
cases. Pairs with size I < contribute to positive ^(ej 1 ) and vice verse. The net E^ determines 
the overall sign of ^(ej 1 ). 

From the dip-recouping shape of $s[Q±(r,cu)] in Fig. [TJ the dip (recouping) region with 
pair size I < I* (I > I*) contributes to positive (negative) ^(ej 1 ), since ^(ej 1 ) cx S(— xt ) (X 
J^ o + drY 2 (r)^s[— Q±(r, to)]/r. The function Y 2 5r oc r(de/dr)5r is proportional to the number 
of pairs present in a ring with radius r and width Sr. Therefore, the sign of ^(e^ 1 ) depends 
on whether there are more pairs with pair size smaller or larger than /* . An interesting point 
to mention is that, upstream vortex motion is not required for negative ^(e^ 1 ); instead, it 
is the asymmetric pair motion in opposite angular directions in Fig. [2] that brings about the 
Hall anomaly. We interpret such result to be that the direction of electric field generated 
by a vortex pair is also pair-size-dependent. This interpretation is illustrated in Fig. HI The 
pair with I < I* polarizes in a such a way that E L points in the direction giving positive 

14 



Hall signal and vice verse. The final sign of ^(e^ 1 ) depends on the net caused by all the 
pairs. 



le-4 




6 8 10 

I —\n(r/a ) 



FIG. 5. (Color online) A plot of Y? versus Z for t = -0.038, 0.021, and 0.030. The vertical arrow 
marks Z*, the position where 9(£7_l) changes sign. 2ttY = 0.1. At t = 0.030, ln(£ + /ao) is indicated 
by the red double-headed arrow. Three representative situations concerning the sign issue of ^(ej 1 ) 
are Y 2 increasing at Z*, Y 2 decreasing at ZJ, and ZJ « ln(£+/ao). 



The weight function Y t 2 , which describes the probability of a pair having pair size r = a^e 1 , 
is plotted in Fig. [5] at different temperatures. The other parameters used are logcJ = —9, 
x = 0.1, and 2nY = 0.1. For t = —0.038, Y? is slightly decreasing at I*, so the dip- 
recouping shape in Sf[(7j_(r, u))] almost cancels each other in the integral in Eq. ( IT7|) . giving 
a weak positive ^(e^ 1 ). If the temperature increases, for example t = 0.021, Y t 2 becomes 
increasing at Z*. More pairs are present in the recouping region and the sign change of 
^(ej 1 ) thus occurs. If the temperature increases further to t — 0.030, the coherence length 
£ + characterizing the maximum pair size becomes comparable to Z* . As a result, only the dip 
part is integrated because the recouping part lies outside of the integration limit in Eq. (|T7|) . 
A strong positive-valued peak then replaces the sign reversal. In summary, temperature 
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controls the pair size distribution, which then determines the sign of ^(e^ 1 ). Since the 
positive slope of at /* only occurs at the high temperature phase, the claim that the sign 
reversal of ^(ej 1 ) is related to KT transition is justified. 

To close this subsection, we discuss the Hall conductivity in the static limit. For low 
temperature phase t < 0, Hall conductivity diverges in zero-frequency limit. Log- log plots of 
^(ej 1 ) and 9ft(— ej 1 ) versus cu using linear fitting show that ^(ej 1 ) oc tu ai and 3?(— ej 1 ) oc cu aR 
with both «r and cq greater than zero but smaller than unity for the temperature range 
t = —0.037 down to —0.15 and frequency range logcJ = —10 to —2. This mean that ej 1 
decreases slower than to when frequency approaches zero, and as a result a± oc ej 1 /{ioS) 
diverges at cu — > 0. 



B. Power dissipation from bound pair contribution 



(a) 12 le-2 t =-0.038 2 le-2 t =0.045 -le-3 <=0.090 
I 

« 0.4 




(h) 14 le-2 logo; = -10 le-2 log^ = -6 le-2 log^ = -2 




0.7 



,0.0, 



J 










-0.1 0.0 0.1 0.2 -0.1 0.0 0.1 0.2 -0.1 0.0 0.1 0.2 

t 



FIG. 6. (Color online) Bound pair contribution. S(— ) as a function of (a) logo; and (b) 



t = T/Tkt — 1- xq = 0.1 and 2ttYq = 0.1. In (a), Sy( — e,, 1 ) increases with frequency when t < 



and is suppressed at small frequency when t > 0. In (b), 9f(— e,, ) has peak structure around Tkt- 
The peak broadens and moves to higher temperature when frequency increases. 
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In this subsection, we discuss the frequency and temperature dependence of ). 
Again, we consider bound pair contribution here. In Fig. EJ plots of ^j(— e^ 1 ) versus (a) 
log uJ and (b) t are presented. The other parameters used are the same as those in Fig. |3j 
In Fig. [6](a), S(— e^ 1 ) increases steadily with frequency at low temperature phase t < 0. 
When t > 0, it is suppressed at small frequency. In Fig. [6](b), 3(— e^ 1 ) first increases 
with temperature, and is suppressed above T KT , resulting in a peak structure. At higher 
frequency, the peak is broadened and shifts to higher temperature. These features can be 
understood under the standard AHNS theory: Above Tkt, the coherence length £ + becomes 
finite, and the integration limits cannot cover the peak structure in O^y) at large I. Thus, 
£y(— e^ 1 ) is suppressed at low frequency for t > in Fig. [6](a), and above some temperatures 
in Fig. [6]^b). Indeed, it is not surprising that the result behaves similarly to its s-wave 
counterpart when we notice that the shape of Q\\ is qualitatively similar to that of g s even 
for finite Xq. 

In zero-frequency limit in the low temperature phase, S(— e^ 1 ) decreases slower than to 
when frequency approaches zero and thus the $t(cr\\) oc £y(— e^ 1 )/^ diverges. This is found 
in log-log plots of e^ 1 ) versus to that S(— e^ 1 ) oc co^ with (3 > but smaller than unity 
using linear fitting from t = —0.037 to —0.15 and logcJ = —10 to —2. However, we cannot 
interpret it to be the divergence of power dissipation because the dissipation also depends 
on the magnitude of the electric field in the superconducting bulk which is supposed to be 
vanishingly small in the static limit. 

C. Total conductivity tensor and free vortex contribution 

Having discussed the result due to bound pair dynamics, we study the total contribution 
from both bound pair and free vortex dynamics in this subsection. In order to discern the 
contribution of bound pair and free vortex dynamics to the total Hall conductivity and 
power dissipation, we plot the total conductivity tensor uLk(T± and uLk^t((T\\) in Figs. [7] 
and M as a function of (a) logcJ and (b) t, as well as the 3(— e^ 1 ), 3?(— ej 1 ), and ^(ej 1 ) due 
purely to the free vortex motion in Fig. [9j We can compare these figures with those which 
take only the bound pair dynamics into account (Figs. [3] and |6]). For the sake of illustration, 
we use C\ = 1 and 6 = 2 for the weight of free vortex contribution. 

For the frequency dependence in Figs. [7] and |8] at temperature t = 0.045, we can see 
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FIG. 7. (Color online) Total ujL^a±_ as a function of (a) logo; and (b) t. The red squares and the 
blue circles represent the real part and imaginary part of u>Lj c a± respectively. C\ = 1 and 6 = 2 
for illustration. The other parameters used are the same as those in Fig. [3l 

free vortex signal emerging at small frequency region. If the temperature increases further 
at t = 0.090, the free vortex signal can merge with and outweigh the bound pair signal 
[compared with Figs. Sandra) t = 0.090]. For t < (not shown), since there is no free 
vortex contribution, we expect the curves are the same as those in low temperature phase 
in Figs. [3] and E(a). As for the temperature dependence in Figs. and E^b), we can see 
for logcJ = —10 the free vortex signal appears at temperature very close to the bound pair 
signal and they almost merge into each other. When the frequency increases (logcU = —6), 
this extra strong signal moves to higher temperature and becomes distinguishable from the 
bound pair signal. The magnitude of the free vortex signal can be seen in Fig. [93 

As discussed in the previous two subsections, 9ft(o"j_), 9?(<7j_), and $t(cr\\) diverge when 
oj — > 0. However, it is only true in low temperature phase. In the case where t > 0, free 
vortex contribution dominates in the static limit 7o/w 3> 1. With such condition, we have 
f k Xf which is purely imaginary. From Eq. (fl9|) . the diagonal part of — Q^xf 1 ) and the 
off-diagonal part of O^x^ 1 ) are given by u/[^ (l + Xq)] and x w/[7o(l + Xq)] respectively, 
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FIG. 8. (Color online) Total uLfJR.(a\\) as a function of (a) logo; and (b) t. The other parameters 
used are the same as those in Fig. [7J 

and therefore both of them are proportional to This shows that a oc xj l /{~ ^ UJ ) is 

purely real, constant in frequency, and diverges when T — > T^ T . As a result, in the high 
temperature phase, both the Hall conductivity and dissipation are protected from divergence 
in the zero-frequency limit but have strong signal when approaching Tkt from above. 

Finally, Fig. [9] shows ^(-e^ 1 ), ^(-ej 1 ), and ^(e^ 1 ) as a function of (a) logo; and (b) 
t including only the free vortex but not the bound pair contribution. The frequency and 
temperature range are chosen to match those in Figs. [7] and El By Eq. ( fT9l) . we see that 
important features in the figure appear around the crossover from 'jo/u ^> 1 to the opposite 
limit 7o/w < 1. In Fig. |9]^a) on the right panel, when observing from small frequency 
region ^q/uj ^> 1, we see that all three quantities increase with frequency. When frequency 
increases further, ^(e^ 1 ) changes sign and then all the quantities decrease in magnitude with 
frequency, which means that it starts to be inefficient for the free vortex to respond to the 
perturbation when frequency is high. In Fig. [9](b), similar features can be observed in the 
temperature domain. Besides, we can see that the main features in the curves emerge at 
a temperature closer to Tkt for small frequency (left panel) than for high frequency (right 
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FIG. 9. (Color online) Free vortex contribution to e |j" 1 )> e J ); an d ^l 6 !*) as a function of 
(a) logcJ and (b) t. The green triangles represent ejT 1 ). The blue circles and the red squares 
represent the negative real part and imaginary part of e^ 1 respectively. The insets show the full 
plot range of ejT ). The other parameters used are the same as those in Fig. [7J 

panel). Although the sign change in ^(e^ 1 ) can also be observed in free vortex picture, it 
does not arise from the response function Q± mentioned before. Instead, it is attributed to 
the pole structure 1 + i^o/u = 0. 

IV. SUMMARY AND REMARK 

In this work, we have generalized AHNS's vortex dynamics in the chiral p-wave supercon- 
ducting state. The behavior of the conductivity tensor near KT transition is investigated. 
We show that the Hall conductivity can be nonzero arising from the vortex-antivortex pair 
unbinding process or from free vortex motion even in the absence of magnetic field. Power 
dissipation is also predicted in this dynamical picture. In low temperature phase, both the 
a± and ^(ctm) diverge in the static limit. But in high temperature phase, the contribution 
from free vortex motion gives finite results in the static limit. We can distinguish the bound 
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pair and free vortex contribution to the total conductivity tensor by comparing the figures 
from different contribution. 



Sign reversal and strong positive peak in u^R,(a±) are illustrated in some suitable frequency 
close to transition, as well as above the transition temperature at fixed frequencies. On the 
other hand, the ojdt(a\\) behaves in a fashion similar to that of the s-wave case. In bound 
pair description, the induced Hall conductivity is strongly influenced by the dip-recouping 
shape of transverse response function which stems from the asymmetric angular 

response of pairs to the driving field. The convective term in the Fokker-Planck equation, 
which originates from the broken T symmetry nature of a chiral p-wave SC, is essential to 
this asymmetry. All these results depend solely on the convective ratio xq, and are valid 
even without applied magnetic field. 

In our work, temperature dependence of the drag coefficients is omitted. Theoretical 
calculation of the drag coefficients has been performed microscopically for SCs^ and for 
and 



SF 3 He (Refs. |3j and l35l ) at temperatures considerably lower than the superconducting 
transition temperature by the use of quasiclassical theory. Generalization of such calcu- 
lation in chiral p-wave pairing state around Tkt will be useful to estimate the strength 
of the convective ratio. Indeed, if we use the formulas of drag coefficients which only 
take into account the Caroli-deGennes-Matricon mode at very low temperature 77W = 
{2iihp°Jm*)uj Q T^ /[{u Q T^Y + 1] and r/' W = (2TThp s /m*)(cu T^) 2 /[(cu t^) 2 + 1]& we can 
arrive at the result = k^Tm* / (2nhp®tuoT^) and = 0, which results in x = 0. How- 
ever, we believe that the effect of temperature and delocalized excitations can contribute to 
the convective ratio. 



Appendix A: Appendix 



In this Appendix, we apply certain approximation to obtain analytic expression of 
S(— e^ 1 ), ^(ej 1 ), and 9R(— ej 1 ) in the bound pair dynamics picture. 
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1. Analytic expression of susceptibility matrix elements 



We first obtain analytic expression for susceptibility matrix elements xt an< ^ Xb by em ~ 
ploying an approximate solution 

, , 1 



9±{r,u) 



1 — iur 2 / (XD) ± ixo 

and taking r(de/dr) to be a constant function of r. We rewrite Eq. (fTTj) in the form 



(Al) 



Xb 



dr— 
if dr 

J a 



9+ + 9- 



dr 



2 

9+ -9 



ix 



g+-g- 



IX 



2 

9+ + 9 



and further approximate the two functions appearing in the above integrands by 



9+ + 9- 
2 

g+-g~ 



1 



Xr 



(r 4 — r) + sgn(xo)r/i5(r — r-y) 



+i [rI 2 S(r — r 2 ) — x^rlidir — r x )] 
rl\5(r — ri) + irl 3 5(r — r 3 ) 
1 



-ix 



Xr 



where 



h = - 



i r , i 

dr- 



1 



6 l (r 4 — r) + sgn(xo)r/i5(r — ri) 



1 



2J r Vl + [ur 2 /(XD) - x ] 2 I + [cur 2 / (\D) + x ]< 
4 + Ax 2 ' 



o 



o 



1 
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7T 

4' 
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dr 



cor 2 /(\D) 
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ur 2 /{\D) 



r V 1 + [wr 2 /(AL>) - x ] 2 1 + [wr 2 /(AL>) + x ] 2 
ur 2 /(\D) ur 2 /{\D) 



2 Jo r V 1 + [wr 2 /(AD) - x ] 2 1 + [wr 2 /(AL>) + x Q \' 



cot 



-1 
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— — 2 + -sgn^xo 
-1 +xi 



7T 
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2' 



(A2) 
(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



In the above expression, r 1; r 2 , r 3 , and r 4 are some characteristic lengths defined as 
follows: ri, r 2 , and r 3 are the peak position of the integrand of Ii, and J 3 respectively. 
r 4 is chosen such that 1/{1 + [wr|/(AD) + x ] 2 } = 1/[2(1 + Xq)]. They are given respectively 

2 V /1 + X 2 +^)V2 / ( V /3 U ;) ) r 2 = r 2 



A-Da/1 + Xq/uj, and 



by r 2 = XD(-1 + x\ 

XD(^/1 + 2x 2 -x )/uj. One more characteristic length worth mentioning is the pair size r* 
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at which £y(£7j_) changes its sign. Within the approximation described above, we find that 
r* = T2 = r^. Then we can obtain the quantity I* = ln(r*/ao) mentioned in the main text. 

Finally, using rde/dr = 4n 3 KoY 2 (r), we obtain analytic expression for the components 
of the susceptibility matrix for £ + > Tj 

xl ~ ?4 - 1 + sgn(x )x vr 4 K Y 2 

(2x 7T \ 

cot" 1 — + ° + -sign x j ir 3 K Yi + i^K Y 2 , (A9) 



(2x 7T \ 

cot_1 1 , 2 + o sign x tt 3 K Y 2 + tx ir 4 K (Y 2 2 - Y 2 ), (A10) 
— 1 + x Z J 

where some shorthand notations are introduced: q = ^(r 4 ) and Yj = Y(rj). From the 
expression above, we notice that x\ (Xb) * s ( an ti-) symmetric in xq as it should be. We can 
also identify the contribution from the dip and recouping part of 3(£?jJ to be Y 2 and Y 2 in 
the imaginary part of Eq. flAlOl) . 



2. Analytic expression of 1 ), ^(e^ 1 ), and e^ 1 ) 

In this subsection, 3(— e^ 1 ), ^(ej 1 ), and 9ft(— ej 1 ) for nonzero Xo are calculated in two 
limiting cases xo <C 1 and xo ^> 1. For ig C 1, we take r| « r| = rf « XD/u, r\ m 
XD/(V3cu), and 

cot" 1 + -sign x « 2x - — ^ + — ^ . (All) 



At the opposite limit xq ^> 1, we take r\ ~ A£)xo(— 1 + \ 2)/u, r\ « rf = r§ « XDxq/uj, 
and 



2xn \ 7T , s , „ 2 



cot ~~ — 2 + o s S n ( x o) ~ Trsgn(xo) + ^-3 - —5 + • ■ • . (A12) 

\ — 1+XqJ 2 X 3Xq 5Xq 

Besides, we neglect the x dependence of ? 4 and Yj in the expansion. Such approximation 
is valid for x <C 1 and all temperature near Tkt because Tj is almost independent of Xq. It 
is also valid for xo ^> 1 in low temperature phase because the renormalization almost go to 
completion at large xo and Yj can be treated as independent of xq oc r 2 . But it is not valid 
for xq ^> 1 in high temperature phase for two reasons: first, Yj is increasing at large Xo and 
second, rj are very likely to be greater than £ + so that the integrations in Eqs. (1A2I) and 



(IA3I) give the trivial results xf, = e (£+) ~~ 1 an d Xb = 0- 
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Expanding using small Xq, the imaginary part of —6m 1 is approximately given by 



« 7fK 7r% 2 (l - 2|x |?I 1 ^ovr 4 F 1 2 ), (A13) 

where we have used 64 ^> K tt a Y£ in the last step. The presence of the second term acts 
as a suppression of dissipation when x® is switched on. If we neglect the correction term 
proportional to \xo\, the expression reduces to the s-wave SF/SCs result. In the opposite 
limit xq 1, the imaginary part of — e^ 1 to the leading order in the smallness of 1/xq is 
given by 



<s -e„ 



(?4 + 2K vr 3 F 1 2 ) 2 + AK 2 ^Yt 
*efK 7r% 2 , (A14) 

assuming €4 3> 2i"f 7r 3 Y 1 2 . Surprisingly, the leading order of 3(— e^ 1 ) at large x has the 
same form as its small Xq counterpart. The only difference is that the characteristic length 



is changed from yXD/u to y XC /(2ttK ui). This can be understood if we recall that the 
shape of O^y) has no qualitative change for finite So- 
Following similar procedure we can obtain analytic expression for ^(ej 1 ) and 9ft(— ej 1 ). 
For small Xq, we obtain ^(e^ 1 ) to the first order in x 



x [?lK o7 r 4 (Y? - Yj) + 4? 4 if 2 7r 7 F 4 4 + Kl^Yj - Kl^YjY?} 

(?1 + K W 4 4 ) 2 

w x 7fK n 4 [(Y? - F 4 2 ) + Ae^K^Y*] . (A15) 

We have used the condition <T 4 S> Kq% a Y 2 in the last line. From the expression, the behavior 
depends on whether Y 2 — Y 2 or 4? 4 ~ 1 i^ 7r 3 F 4 4 dominates. From numerical results, Y 2 — 
Y 4 2 dominates at a wide range of temperature in the low temperature phase. Only when 
approaching the transition temperature from below and away from the vicinity w^ 2 ./ {\D) ~ 
1, the Y£ becomes more significant. On the other hand, 

2x (? A K ^Yi +? 4 ^ W 4 4 - Kl-K n Y! -e 4 Kfr*YiY?) 



(ei + Kl^Yif 

2x efK ir 3 Yl (A16) 
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We have used 3> Kqtc'Y'^ in the second line. For large xo, we assume 3> AKq^Y^ and 
obtain to the leading order 



sgn(xo) 



? 4 + 2K 7r 3 F 1 2 )(?2 + 4?4^o7r 3 ^ 2 + 4Zf 2 7r 6 Y 1 4 + 4K 2 vr 8 F 1 4 ) 
« sgn(x )2K 2 7r 8 ?7 3 >i\ (A17) 

and 

~ Sgn^oj^ + ^ K ^Zy2 + 4 ^2 7r 6 F 4 + AK 2^ 

« sgn^o)?! 2 ^^ 4 ^ 2 . (A18) 

It has to be mentioned that, although we have used the condition q 3> Kq-k Y? extensively, 
it is only valid in the temperature range t < 0.035. Above such temperature, numerical 
results show that Kqtt^Y? can be greater than q at small frequency region. 
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